***********阿特金森指数******************
sum 总收入
replace 税前指数1=总收入/106647.8 
replace 税前指数1=税前指数1^(0.5)  //系数为0.5 
egen sum_税前指数1=sum(税前指数1)
replace 税前指数1=sum_税前指数1/8727
replace 税前指数1=税前指数1^2
replace 税前指数1=1-税前指数1

replace 税前指数2=log(总收入) //系数为1
egen sum_税前指数2=sum(税前指数2)
replace 税前指数2=sum_税前指数2/8727
replace 税前指数2=2.71828^(税前指数2)  //ye
replace 税前指数2=1-(税前指数2/106647.8 )

replace 税前指数3=总收入/106647.8 
replace 税前指数3=税前指数3^(-1)  //系数为2
egen sum_税前指数3=sum(税前指数3)
replace 税前指数3=sum_税前指数3/8727
replace 税前指数3=税前指数3^(-1) 
replace 税前指数3=1-税前指数3

replace 税后总指数1=总收入-总税额  //总税额
sum 税后总指数1
replace 税后总指数1=税后总指数1/101464.4 
replace 税后总指数1=税后总指数1^(0.5)  //系数为0.5 
egen sum_税后总指数1=sum(税后总指数1)
replace 税后总指数1=sum_税后总指数1/8727
replace 税后总指数1=税后总指数1^2
replace 税后总指数1=1-税后总指数1

replace 税后总指数2=总收入-总税额  //总税额
replace 税后总指数2=log(税后总指数2) //系数为1
egen sum_税后总指数2=sum(税后总指数2)
replace 税后总指数2=sum_税后总指数2/8727
replace 税后总指数2=2.71828^(税后总指数2)  //ye
replace 税后总指数2=1-(税后总指数2/101464.4  )


replace 税后总指数3=总收入-总税额  //总税额
replace 税后总指数3=税后总指数3/101464.4 
replace 税后总指数3=税后总指数3^(-1)  //系数为2
egen sum_税后总指数3=sum(税后总指数3)
replace 税后总指数3=sum_税后总指数3/8727
replace 税后总指数3=税后总指数3^(-1)
replace 税后总指数3=1-税后总指数3

replace 税后房贷指数1=总收入-房贷价值  //房贷价值
replace 税后房贷指数1=税后房贷指数1/104964.9 
replace 税后房贷指数1=税后房贷指数1^(0.5)  //系数为0.5 
egen sum_税后房贷指数1=sum(税后房贷指数1)
replace 税后房贷指数1=sum_税后房贷指数1/8727
replace 税后房贷指数1=税后房贷指数1^2
replace 税后房贷指数1=1-税后房贷指数1

replace 税后房贷指数2=总收入-房贷价值  //房贷价值
replace 税后房贷指数2=log(税后房贷指数2) //系数为1
egen sum_税后房贷指数2=sum(税后房贷指数2)
replace 税后房贷指数2=sum_税后房贷指数2/8727
replace 税后房贷指数2=2.71828^(税后房贷指数2)  //ye
replace 税后房贷指数2=1-(税后房贷指数2/104964.9)

replace 税后房贷指数3=总收入-房贷价值  //房贷价值
replace 税后房贷指数3=税后房贷指数3/104964.9
replace 税后房贷指数3=税后房贷指数3^(-1)  //系数为2
egen sum_税后房贷指数3=sum(税后房贷指数3)
replace 税后房贷指数3=sum_税后房贷指数3/8727
replace 税后房贷指数3=税后房贷指数3^(-1)
replace 税后房贷指数3=1-税后房贷指数3

replace 税后房贷福利指数1=总收入-房贷价值+家庭房贷价值福利  // 家庭房贷价值福利
replace 税后房贷福利指数1=税后房贷福利指数1/106647.8
replace 税后房贷福利指数1=税后房贷福利指数1^(0.5)  //系数为0.5 
egen sum_税后房贷福利指数1=sum(税后房贷福利指数1)
replace 税后房贷福利指数1=sum_税后房贷福利指数1/8727
replace 税后房贷福利指数1=税后房贷福利指数1^2
replace 税后房贷福利指数1=1-税后房贷福利指数1

replace 税后房贷福利指数2=总收入-房贷价值+家庭房贷价值福利  //家庭房贷价值福利
replace 税后房贷福利指数2=log(税后房贷福利指数2) //系数为1
egen sum_税后房贷福利指数2=sum(税后房贷福利指数2)
replace 税后房贷福利指数2=sum_税后房贷福利指数2/8727
replace 税后房贷福利指数2=2.71828^(税后房贷福利指数2)  //ye
replace 税后房贷福利指数2=1-(税后房贷福利指数2/106647.8)


replace 税后房贷福利指数3=总收入-房贷价值+家庭房贷价值福利  // 家庭房贷价值福利
replace 税后房贷福利指数3=税后房贷福利指数3/106647.8
replace 税后房贷福利指数3=税后房贷福利指数3^(-1)  //系数为2
egen sum_税后房贷福利指数3=sum(税后房贷福利指数3)
replace 税后房贷福利指数3=sum_税后房贷福利指数3/8727
replace 税后房贷福利指数3=税后房贷福利指数3^(-1)
replace 税后房贷福利指数3=1-税后房贷福利指数3

replace 税后未成年福利指数1=总收入-房贷价值+家庭未成年房贷价值福利  // 家庭未成年房贷价值福利
replace 税后未成年福利指数1=税后未成年福利指数1/106647.8
replace 税后未成年福利指数1=税后未成年福利指数1^(0.5)  //系数为0.5 
egen sum_税后未成年福利指数1=sum(税后未成年福利指数1)
replace 税后未成年福利指数1=sum_税后未成年福利指数1/8727
replace 税后未成年福利指数1=税后未成年福利指数1^2
replace 税后未成年福利指数1=1-税后未成年福利指数1

replace 税后未成年福利指数2=总收入-房贷价值+家庭未成年房贷价值福利  //家庭未成年房贷价值福利
replace 税后未成年福利指数2=log(税后未成年福利指数2) //系数为1
egen sum_税后未成年福利指数2=sum(税后未成年福利指数2)
replace 税后未成年福利指数2=sum_税后未成年福利指数2/8727
replace 税后未成年福利指数2=2.71828^(税后未成年福利指数2)  //ye
replace 税后未成年福利指数2=1-(税后未成年福利指数2/106647.8)

replace 税后未成年福利指数3=总收入-房贷价值+家庭未成年房贷价值福利  // 家庭未成年房贷价值福利
replace 税后未成年福利指数3=税后未成年福利指数3/106647.8
replace 税后未成年福利指数3=税后未成年福利指数3^(-1)  //系数为2
egen sum_税后未成年福利指数3=sum(税后未成年福利指数3)
replace 税后未成年福利指数3=sum_税后未成年福利指数3/8727
replace 税后未成年福利指数3=税后未成年福利指数3^(-1)
replace 税后未成年福利指数3=1-税后未成年福利指数3

gen 社会福利效应总指数1=((101464.4*(1-税后总指数1))/(106647.8*(1-税前指数1)))^0.5  //系数为0.5
gen 社会福利效应总指数2=(101464.4*(1-税后总指数2))/(106647.8*(1-税前指数2))  //系数为1
gen 社会福利效应总指数3=((101464.4*(1-税后总指数3))/(106647.8*(1-税前指数3)))^(-1)  //系数为2

gen 社会福利房贷价值指数1=((104964.9*(1-税后房贷指数1))/(106647.8*(1-税前指数1)))^0.5  //系数为0.5
gen 社会福利房贷价值指数2=(104964.9*(1-税后房贷指数2))/(106647.8*(1-税前指数2))  //系数为1
gen 社会福利房贷价值指数3=((104964.9*(1-税后房贷指数3))/(106647.8*(1-税前指数3)))^(-1)  //系数为2

gen 社会福利房贷福利指数1=((106647.8*(1-税后房贷福利指数1))/(106647.8*(1-税前指数1)))^0.5  //系数为0.5
gen 社会福利房贷福利指数2=(106647.8*(1-税后房贷福利指数2))/(106647.8*(1-税前指数2))  //系数为1
gen 社会福利房贷福利指数3=((106647.8*(1-税后房贷福利指数3))/(106647.8*(1-税前指数3)))^(-1)  //系数为2

gen 社会福利未成年福利指数1=((106647.8*(1-税后未成年福利指数1))/(106647.8*(1-税前指数1)))^0.5  //系数为0.5
gen 社会福利未成年福利指数2=(106647.8*(1-税后未成年福利指数2))/(106647.8*(1-税前指数2))  //系数为1
gen 社会福利未成年福利指数3=((106647.8*(1-税后未成年福利指数3))/(106647.8*(1-税前指数3)))^(-1)  //系数为2

sum 总收入 if 收入阶级==2
gen 税前指数1=总收入/109101.6    if 收入阶级==2
replace 税前指数1=税前指数1^(0.5) if 收入阶级==2 //系数为0.5 
egen sum_税前指数1=sum(税前指数1) if 收入阶级==2
replace 税前指数1=sum_税前指数1/2133  if 收入阶级==2
replace 税前指数1=税前指数1^2 if 收入阶级==2
replace 税前指数1=1-税前指数1 if 收入阶级==2

gen 税前指数2=log(总收入)  if 收入阶级==2 //系数为1
egen sum_税前指数2=sum(税前指数2) if 收入阶级==2
replace 税前指数2=sum_税前指数2/2133  if 收入阶级==2
replace 税前指数2=2.71828^(税前指数2) if 收入阶级==2 //ye
replace 税前指数2=1-(税前指数2/109101.6 ) if 收入阶级==2

gen 税前指数3=总收入/109101.6   if 收入阶级==2
replace 税前指数3=税前指数3^(-1) if 收入阶级==2 //系数为2
egen sum_税前指数3=sum(税前指数3) if 收入阶级==2
replace 税前指数3=sum_税前指数3/2133   if 收入阶级==2
replace 税前指数3=税前指数3^(-1) if 收入阶级==2
replace 税前指数3=1-税前指数3 if 收入阶级==2


gen 税后房贷福利指数1=总收入-房贷价值+家庭房贷价值福利  if 收入阶级==4 // 家庭未成年房贷价值福利
sum 税后房贷福利指数1 if 收入阶级==4 
replace 税后房贷福利指数1=税后房贷福利指数1/337496.5    if 收入阶级==4
replace 税后房贷福利指数1=税后房贷福利指数1^(0.5) if 收入阶级==4 //系数为0.5 
egen sum_税后房贷福利指数1=sum(税后房贷福利指数1) if 收入阶级==4
replace 税后房贷福利指数1=sum_税后房贷福利指数1/669 if 收入阶级==4
replace 税后房贷福利指数1=税后房贷福利指数1^2 if 收入阶级==4
replace 税后房贷福利指数1=1-税后房贷福利指数1 if 收入阶级==4

gen 税后房贷福利指数2=总收入-房贷价值+家庭房贷价值福利 if 收入阶级==4 //家庭未成年房贷价值福利
replace 税后房贷福利指数2=log(税后房贷福利指数2) if 收入阶级==4 //系数为1
egen sum_税后房贷福利指数2=sum(税后房贷福利指数2) if 收入阶级==4
replace 税后房贷福利指数2=sum_税后房贷福利指数2/669 if 收入阶级==4
replace 税后房贷福利指数2=2.71828^(税后房贷福利指数2) if 收入阶级==4 //ye
replace 税后房贷福利指数2=1-(税后房贷福利指数2/337496.5) if 收入阶级==4

gen 税后房贷福利指数3=总收入-房贷价值+家庭房贷价值福利 if 收入阶级==4 // 家庭未成年房贷价值福利
replace 税后房贷福利指数3=税后房贷福利指数3/337496.5    if 收入阶级==4
replace 税后房贷福利指数3=税后房贷福利指数3^(-1) if 收入阶级==4 //系数为2
egen sum_税后房贷福利指数3=sum(税后房贷福利指数3) if 收入阶级==4
replace 税后房贷福利指数3=sum_税后房贷福利指数3/669 if 收入阶级==4
replace 税后房贷福利指数3=税后房贷福利指数3^(-1) if 收入阶级==4
replace 税后房贷福利指数3=1-税后房贷福利指数3 if 收入阶级==4

gen 社会福利房贷福利指数1=((109897.3 *(1-0.0024))/(109101.6*(1-0.0022)))^0.5  if 收入阶级==2 //系数为0.5
gen 社会福利房贷福利指数2=(109897.3 *(1-0.0049))/(109101.6*(1-0.0044))  if 收入阶级==2 //系数为1
gen 社会福利房贷福利指数3=((109897.3  *(1-0.0097))/(109101.6*(1-0.0088)))^(-1)  if 收入阶级==2 //系数为2

*******无分配**********************
sum 总收入 

sum 自有套数 

sum number 

sum 人均建筑面积 

sum 总房贷


sum 房贷价值

***********公共服务均等分配*************
egen sum_房贷价值=sum(房贷价值)
egen sum_number=sum(number)
gen 人均房贷价值福利=sum_房贷价值/sum_number
gen 家庭房贷价值福利=number*人均房贷价值福利
gen 房贷价值净福利=家庭房贷价值福利-房贷价值
sum 房贷价值净福利 if 房贷价值净福利>0
sum 房贷价值净福利 if 房贷价值净福利<=0
sum 总收入 if 房贷价值净福利>0
sum 总收入 if 房贷价值净福利<=0
sum 自有套数 if 房贷价值净福利>0
sum 自有套数 if 房贷价值净福利<=0
sum number if 房贷价值净福利>0
sum number if 房贷价值净福利<=0
sum 人均建筑面积 if 房贷价值净福利>0
sum 人均建筑面积 if 房贷价值净福利<=0
sum 总房贷 if 房贷价值净福利>0
sum 总房贷 if 房贷价值净福利<=0
egen sum_房贷价值净福利=sum(房贷价值净福利)
***********仅分配给有未成年子女的家庭*************
egen sum_房贷价值=sum(房贷价值)
egen sum_未成年=sum(未成年)
gen 人均未成年房贷价值福利=sum_房贷价值/sum_未成年
gen 家庭未成年房贷价值福利=人均未成年房贷价值福利*未成年
gen 未成年房贷价值净福利=家庭未成年房贷价值福利-房贷价值
sum 未成年房贷价值净福利 if 未成年房贷价值净福利>0
sum 未成年房贷价值净福利 if 未成年房贷价值净福利<=0
sum 总收入 if 未成年房贷价值净福利>0
sum 总收入 if 未成年房贷价值净福利<=0
sum 自有套数 if 未成年房贷价值净福利>0
sum 自有套数 if 未成年房贷价值净福利<=0
sum number if 未成年房贷价值净福利>0
sum number if 未成年房贷价值净福利<=0
sum 人均建筑面积 if 未成年房贷价值净福利>0
sum 人均建筑面积 if 未成年房贷价值净福利<=0
sum 总房贷 if 未成年房贷价值净福利>0
sum 总房贷 if 未成年房贷价值净福利<=0
egen sum_未成年房贷价值净福利=sum(未成年房贷价值净福利)
************净福利决定因素********************
gen 总面积=一套面积+二套面积+三套面积+四套面积+五套面积+六套面积
replace 是否拆迁=0 if 是否拆迁==2
replace 住房=1 if j1001==1
sum  房贷价值净福利 number 总收入 总现值 自有套数 人均建筑面积  总储蓄 住房 是否拆迁 房价 税率
gen ln房贷价值净福利=ln(房贷价值净福利)
gen ln总收入=ln(总收入)
gen ln总现值=ln(总现值)
gen ln总储蓄=ln(总储蓄)
gen ln房价=ln(房价)
gen ln人均建筑面积=ln(人均建筑面积)
gen 无减免福利=-房贷价值
sum 无减免福利
sum 房贷价值净福利
sum 未成年房贷价值净福利
sum 税率
reg 无减免福利 number 自有套数 总面积 总现值 总收入 总储蓄  基础消费  住房 是否拆迁 税率, r 
reg 房贷价值净福利 number 自有套数 总面积 总现值 总收入 总储蓄  基础消费  住房 是否拆迁 税率, r 
reg 未成年房贷价值净福利 number 自有套数 总面积 总现值 总收入 总储蓄 基础消费 住房 是否拆迁 税率, r

preserve
keep if 省份1==31 | 省份1==32 | 省份1==33 | 省份1==35 | 省份1==37 | 省份1==11 | 省份1==12 | 省份1==13 | 省份1==44 | 省份1==46 | 省份1==21 | 省份1==22 | 省份1==23
reg 未成年房贷价值净福利 number 自有套数 总面积 总现值 总收入 总储蓄 基础消费 住房 是否拆迁 税率, r
restore 

preserve
keep if 省份1==14 | 省份1==34 | 省份1==36 | 省份1==41 | 省份1==42 | 省份1==43 | 省份1==45  
reg 未成年房贷价值净福利 number 自有套数 总面积 总现值 总收入 总储蓄 基础消费 住房 是否拆迁 税率, r
restore 

preserve
keep if 省份1==61 | 省份1==62 | 省份1==63 | 省份1==64 | 省份1==50 | 省份1==51 | 省份1==52 | 省份1==53
reg 未成年房贷价值净福利 number 自有套数 总面积 总现值 总收入 总储蓄 基础消费 住房 是否拆迁 税率, r
restore 


tab 房价
preserve
keep if 房价>7000
reg 未成年房贷价值净福利 number 自有套数 总面积 总现值 总收入 总储蓄 基础消费 住房 是否拆迁 税率, r
restore 

preserve
keep if 房价<=7000 
reg 未成年房贷价值净福利 number 自有套数 总面积 总现值 总收入 总储蓄 基础消费 住房 是否拆迁 税率, r
restore 

